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Introduction 

The  work  conducted  under  this^  contract  constitutes  a  data-processing 
and  data  analysis  part  of  the  AFGL  SCRIBE  (Stratos])heric  Cryogenic  Inter¬ 
ferometer  Balloon  Experiment J  program.  Our  participation  may  be  summarized 
into  three  major  categories: 

(1)  extraction  of  the  interferogram  data  from  the  PCM  telemetry 
data  tape; 

(2)  processing  of  the  interferogram  data  for  spectral  retrieval; 
and 

(3)  analysis  of  the  obtained  atmospheric  emission  data. 

The  first  successful  balloon  flight  took  place  during  the  month  of 
October  1980.  Our  effort  prior  to  the  flight  for  preparing  various 
electronic  hardware,  as  well  as  necessary  computer  software  for  the  data 
processing,  was  proved  to  be  effective  when  they  were  applied  for 
processing  the  flight  data  which  became  available.  We  were  able  to  extract 
the  data  without  encountering  difficulties  and  to  analyze  them  successfully 

Our  scientific  report  entitled  "SCRIBE  I  Data  Analysis"  describes  the 
technical  aspect  of  our  effort.^  A  paper  entitled  "Measurement  of  atmo¬ 
spheric  emission  using  a  balloon-borne  cryogenic  Fourier  spectrometer"  was 
published  in  the  Proceedings  of  the  1981  International  Conference  on 
Fourier  Transform  Infrared  Spectroscopy.^ 

^H.  Sakai,  "SCRIBE  I  Data  Analysis,"  AFGL-TI<-bl-0:i29,  April  198i  • 

^A  copy  is  supplied  with  this  final  report. 
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The  scheme  developed  under  this  contract  proved  its  effectiveness  for 
the  data  processing  performed  on  the  SCRIBE  II  data  of  October  I98I.  We  are 


able  to  process  more  than  10  interferogram  data  a  day. 

In  this  final  report,  we  focus  our  attention  to  the  operational  aspect 
of  our  processing  scheme. 


SCRIBE  Program 

A  main  goal  of  the  SCRIBE  program  is  to  obtain  the  atmospheric  emission 
spectral  data  at  a  balloon  altitude  of  30,000  m.  The  spectrometer's  line 
of  sight  is  set  at  a  direction  within  a  few  degrees  above  the  horizon.  The 
spectrometer  used  for  observation  is  a  cat's  eye  retroflector  Michelson 
interferometer  operated  at  liquid  nitrogen  temperature  with  a  Ge:Hg  or 
GE;Cu  detector.  The  infrared  radiation  which  we  want  to  observe  is  the 
various  molecular  bands  of  the  atmospheric  species  emitted  under  the  local 
thermal  equilibrium  condition  at  230®K,  which  is  the  ambient  temperature  at 
30  Km  height. 

The  molecular  density  of  major  molecular  species  at  the  altitude  is 
given  in  Table  I. 


Table  I 


n 

o 
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Molecular  Density  per  cm 
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at 


30  Km  Altitude 
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Accepting  the  atmospheric  scale  height  of  8  Km,  we  can  expect  the 


total  molecular  column  density  N  per  cm  approximate] y  given  by 


N  =  n^(8  X  lO^)/cos0, 

where  9  is  the  elevation  angle  of  the  sight  line.  At  the  altitude  of  30  Km, 
the  ambient  pressure  is  an  order  of  10  mb.  Botli  the  collisJonal  broadening 
process  and  the  Doppler  broadening  process  are  equally  important  at  the 
altitude  of  30  Km.  The  molecular  line  profile  k(o)  with  the  Gaussian  form 
is  given  by 


k(a) 


-  In^i - 

“d 


where  S  is  the  total  line  strength  which  is  computed  by  multiplying  the 
total  colimn  density  to  the  absorption  strength  G^, 

G  =  s  N  =  /k(o)do. 
o 

The  Doppler  line  width  is  a  function  solely  determined  by  the  temperature 
T,  the  molecular  mass  m,  and  the  line  center  frequency  0^: 

a  /3M 

d  c  ^  m  ' 

The  emission  radiance  of  a  molecular  line  calculated  under  the  single  layer 
assumption  is  given  by 

E(a)  =  B(o,T)  [1  -  exp(-  e  )], 

d 

where  B{o,T)  is  the  blackbody  radiation  function  at  a  and  T. 

If  the  molecular  line  broadening  is  controlled  by  the  collisional 
process,  the  profile  is  Lorentzian: 
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k(o)  = 


Sa 

c _ 

r/  s2  ^  2,  ' 

ir[(a-a  )  +  a  J 

o  c 

Again  we  get 

S  =  /k(o)da. 

The  emission  radiance  of  a  molecular  line  calculated  under  the  same 
assumption  is  given  by 

S  Ita 

E(a)  =  B(a,T)  [l  -  exp  -  ( - - 5}]. 

The  power  available  for  observation  of  the  molecular  line  is  then 

determined  by  the  radiance  power  E  times  the  interferometer  energy  collecting 

-U  2 

power  which  is  approximately  1.2  x  10  steradian  cm  for  our  case.  Taking 
the  Doppler  width  of  O.OOO8  cm  ^  and  the  blackbody  radiance  of  5  x  10  ^  W/ 

2  ““j. 

steradian  cm  cm  ,  the  energy  power  for  observation  of  a  single  molecular 

-1  -12 
line  at  the  700  cm  region  is  found  to  be  an  order  of  10  W,  imposing  two 

basic  requirements  for  the  measurement  sensitivity.  The  detector  must  be 

extremely  sensitive  and  the  interferometer  itself  must  be  cooled  to  liquid 

nitrogen  temperature. 

The  interferometer  is  driven  at  a  rate  of  30  sec  to  make  a  single 
complete  scanning  cycle.  The  inter ferogreun  data  acquisition  period  lasts 
approximately  25  sec.  The  interferogram  data  are  transmitted  through  the 
PCM  telemetry  link  to  the  ground  station  at  a  rate  of  11,000  T2-bit  words 
per  second.  The  received  telemetry  signal  is  recorded  on  an  analog  tape 
fPCM  tape]  running  at  60  ips.  Together  with  the  interferogram  signal  in  a 
digital  sequence  are  the  PCM  synchronization  code,  the  other  housekeeping 
signals,  the  time  code,  etc.  The  time  code  provides  the  information  on  the 
time  of  recording  with  the  Greenwich  standard  time. 
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Extraction  of  the  interferoRram  data  from  the  PCM  tape 
The  operation  is  performed  using  an  analog  TAPE  recorder /playback 
machine,  a  time  code  generator,  a  PCM  signal  decoder,  a  home-made  interface 
electronics,  a  PDPll  computer  with  a  RKll  disk  and  a  TMll  magnetic  tape 
transport.  The  operational  sequence  is  detailed  below.  The  PDPll  computer 
is  operated  under  the  RTll  system  software.  The  software  written  for  the 
data  transfer  between  the  DRUB  parallel  input  to  the  RKll  disk  is  listed 
in  our  Scientific  Report  I. 

(1)  Mount  the  PCM  tape  on  a  playback  unit.  Connect  the  PCM  signal  output 
to  a  PCM  decoder  and  the  time-code  output  to  a  time-code  generator. 

(2)  Play  back  the  PCM  tape.  Set  the  PCM  decoder  for  a  proper  synch  word 
and  a  proper  frame  word.  Observe  the  bit  pattern  displayed  on  the 
PCM  decoder  unit.  The  time-code  generator  displays  the  time  recorded 
on  the  PCM  tape. 


(3)  Start  the  PDPll.  Bring  the  keyboard  monitor  up.  Set  the  switch 

register  to  105-  Mount  a  magnetic  tape  to  the  TMll  tape  transport. 
Initialize  the  magnetic  tape  by  typing  MT:/Z  under  PIP.  Bring  back 
KB  monitor. 

(ii)  Place  the  PCM  tape  at  a  proper  position  by  checking  the  time  code. 
Start  the  PCM  tape. 

(5)  Wait  until  the  signal  indicates  the  retrace  mode  (bit  pattern  OOOOOO). 

on  the  keyboard.  The  computer  responds  a  bell  signal 


Type 


R  kj 


indicating  the  program  is  ready.  If  the  bell  signal  comes  after 
completion  of  the  retrace  mode,  go  back  to  KB  monitor  and  repeat  the 
sequence  (I4)  and  (5)« 


(6)  If  everything  is  OK,  the  computer  signals  the  bell  sound  after  it 


transfers  an  interferogram  data  extracted  from  the  PGM  tape.  During 
this  period,  no  synch  loss  should  occur.  After  the  bell  sound,  the 
computer  prints  "DONE."  Stop  the  PCM  tape. 

(7)  Bring  back  KB  monitor.  Run  PIP  by  typing  R  PIP^.  After  USR  responds 
by  printing  *,  type  Mr0:LFN.DAT<RKl:PCMRKl.DAT/M=l^.  Observe  the  data 
transfer  from  RKl : PCMRKl . DAT  to  MT0:LFN.DAT  (LFN;  identification  file 
name) . 

(8)  Repeat  (k)  through  (T)  until  the  MT  tape  is  full.  [Normally  7  data 
transfers  are  required. ] 

(9)  Demount  the  MT  tape  from  the  TMll  tape  transport.  Put  a  g\mimed  label 
on  the  MT  tape  for  its  identification.  The  current  label  is  shown 
below.  It  is  preferred  to  have  another  gummed  label,  indicating  all 


Because  of  a  convenience  provided  by  the  CDC  computer  system  we  have 
at  the  University  of  Massachusetts,  in  particular  because  of  its  convenient 
interactive  featiire,  both  the  spectral  recovery  and  the  spectral  analysis 
are  performed  using  our  central-site  computing  facility. 

The  input  data  for  the  spectral  recovery  process  is  the  raw  interfero- 
gram  data  recorded  on  a  7-traok  magnetic  tape  under  the  control  of  the  RTll 
system  software.*  The  logical  flow  was  described  in  our  scientific  report. 

« 

The  tape  has  a  7-track,  odd-parity  format  recorded  with  800  BPI  density.  A 
single  data  information  consists  of  three  files  separated  by  the  EOF  mark, 
the  directory,  the  data,  and  a  single  empty  file  ended  by  the  EOF. 
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The  decoding  software  was  slightly  modified  from  the  source  listing  in  our 
scientific  report  to  the  one  in  the  Appendix  for  improvement  of  its  operating 
efficiency.  The  phase  correction,  the  Fourier  transformation  and  the  CRT 
plot  software  were  those  developed  by  us  for  other  purposes,  and  their 
documentation  would  appear  elsewhere.  The  processing  scheme  is  sequenced 
as  follows: 

(1)  Run  Program  A  shown  in  Figure  1  with  a  data  tape  ARFII2.  The  PCMDCC 
would  have  a  file  structure  shown  in  Figure  2. 

(2)  Use  COPYBF  utility  program  to  copy  the  data  to  TAPF9. 

(3)  Run  the  program  on  TREAD.  The  output  result  is  on  TAPE2. 

(U)  Move  the  data  on  TAPE2  to  TAPES. 

Return,  TAPE2,  TAPE6,  and  TAPE9. 

( 5 )  Rewind  TAPES . 

Run  the  program  on  SPCXX. 

The  sequences  (2)  through  (5)  are  shown  in  Figure  3. 

(6)  Inspect  the  printed  numbers.  An  example  of  the  printout  is  shown  in 
Figure  U. 

Find  where  the  gain  switch  occurs,  and  find  M  and  N. 

Type  M  and  N  in  (215)  FORMAT 
The  output  is  on  TAPEl. 

(7)  Rewind  TAPEl. 

Run  LFILFIX. 

Type  all  control  parameters  as  asked. 

(8)  Return,  TAPEl 
Rewind,  TAPE2 

Define  TAPES  for  a  direct  access  permanent  file. 


Run  Li,  L2,  L3,  and  LU  consecutively. 

The  output  on  TAPE3  is  a  spectrum  recorded  in  (1615)  format.  It  is 
normalized  to  1000  at  the  highest  peak.  See  the  example  shown  in 
Figtire  k  for  the  sequence  (7)  and  (8). 

(9)  Go  back  to  (l ) . 
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Program  A 


WI45A,  (JOB  HASH  code) 

ACCOUNT  . 

VSN,  TAPEl  =  ARFTI2. 

LABEL  (TAPEl.,  MT,  D=800,  LB±KU,  P*S,  PO=R) 

Define,  PCMDCC. 

COPYBF,  TAPEl,  PCMDCC,  n  /  n=3N-l  when?  N  is  Itie  number  of 

interferograin  dat.i 

DAYFILE,  L=SKDFILE. 

Replace,  SKDFILE. 

EXIT. 

DAYFILE,  L=SKDFILE. 

Replace,  SKDFILE. 

/EOR/ 

/EOI/ 


Figure  1 
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Directory 

EOF 

Data  I 

EOF 

EOF 

Directory 

1 

2 

F  Data  II 


/OEI » TREAD fSKXX 
/FTN»I=aREADi-L-0 

.271  CP  SECONDS  COMPILATION  TIME 
/FTN>  J>SKXX»L==OrD=LQ 

.172  CP  SECONDS  COMPILATION  TIME 
/GETfLEVN^LFlLFlX 
/COPYDRrLEyNrLl 
COPY  COMPLETE. 

/C0PYBRrLEVNrL2r6 
COPY  COMPLETE. 

/COPYDRrLE  VNrLar  7 
COPY  COMPLETE. 

/C0PYDR»LEVNrL4r5 
COPY  COMPLETE. 

COPYBF.PCMDCC.X.n. 

COPYBF , PCMDCC , TAPE9. 

Rewind,TAPE9. 

LCO 

Rewind, TAPE2, 

COPYCR ,TAPE2 , TAPES 

Ret  urn , TAPE? , TAPES , TAPE9 . 

Rewind, TAPES, 


Figure  3 


Spectral  Analysis 

The  SCRIBE  program  is  a  technical  challenge,  attempting  to  achieve  a 
good  measiirement  on  the  atmospheric  emission  data  which  hitherto  have  been 
known  to  us  rather  poorly.  Based  on  the  primary  goal  set  for  the  program, 
our  analysis  is  designed  to  spot  shortcomings  which  may  exist  on  either 
the  observed  data  or  the  presently  existing  theoretical  model,  without 
requiring  an  elaborate  computational  effort.  In  other  words,  our  effort  is 
focused  to  survey  the  observed  data  using  a  theoretically  computed  model 
which  would  provide  a  radiance  figiire  within  ^  20%  accuracy.  The  consider¬ 
ation  dictated  to  us  for  our  analysis  effort  is  a  cost  and  a  speed  of 
computation.  The  balloon  data  would  be  obtained  anyvay  in  a  rather 
difficult  circumstance  with  parameters  which  are  determined  with  limited 
accuracies. 

Upon  considering  these  points,  we  decided  to  synthesize  the  theoretical 
spectrum  primarily  with  an  emphasis  on  the  computational  speed.  At  the  same 
time  the  computational  accxxracy  suffered  to  some  extent,  in  order  to  achieve 
improved  computational  efficiency.  Nonetheless,  the  spectral  result 
obtained  in  the  computation  would  serve  to  provide  a  reference  for  survey 
effort.  We  designed  the  computational  algorithm  to  achieve  a  speed  suitable 
foi-  interactive  study. 

Tiio  atmospheric  line  1 isting  compiled  by  the  AFGL  forms  a  theoretical 
base  for  our  computation.  From  the  beginning,  we  did  not  attempt  to 
generate  an  accurate  line  profile  in  oui-  computation,  because  both  factors, 
insufficient  .'.pectrai  resolution  and  blending  of  lines,  would  conceal 
individual  line  profiles  in  the  observed  data.  We  did  not  attempt  to  use 


IP 


the  AFGL  FASCOD^  for  our  computation  because  of  the  reason  cited  above. 

The  molecular  absorptance  A(o)  at  wavenumber  o  is  calculated  by 

A(o)  =  1  -  exp  -  {Zk  (o)}, 

n 

n 

where  is  the  line  absorbance  of  a  line  n.  The  optical  depth  of  the 

molecular  line  with  which  we  deal  in  our  caiculation  is  in  general  relatively 
thin.  The  spectral  resolution  in  the  observed  data  is  insufficient  for 
allowing  us  to  observe  the  individual  line  profile  without  a  large- 
distortion,  Under  these  conditions,  the  absorptance  averaged  over  a 
spectral  resolution  much  wider  than  the  individual  line  width,  can  be 
calculated  by 


=  1  -  exp  -  <{5;k^(o))>^^. 
n 

Tile  averaging  operation  is  taken  over  t)ie  resolution  widtii;  i.e.,  the 
absorbance  Zk^(o)  is  averaged  over  the  lines  which  fall  within  the  resolu¬ 
tion  width  Ao.  If  n  lines  are  in  Aa,  the  average  absorbance  <Ik^(o)>^y  is 

n 


computed  by 


n  n 

where  is  the  line  strength  of  the  line  n. 

The  energy  is  absorbed  by  the  molecular  system  from  the  radiation 
field  at  a  rate  specified  by  the  absorptance  along  the  absorption  path. 
The  absorbed  energy  by  the  molecular  system  is  redistributed  within  the 
entire  system.  If  a  thermal  equilibrium  is  established  at  temperature  T 
the  molecular  system  and  the  radiation  field  are  in  thermal  equilibi ium 
at  that  temperature.  If  the  molecular  system  is  capable  of  absorbing 


^G.A.  Clough,  F.X.  Kneizis,  and  J.H.  Chetwind,  AFGL-TR-77-016U . 
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the  radiation  energy  by  the  absorptance  A(a),  the  same  molecular  system 
establishes  a  thermal  equilibrium  with  the  radiation  field  which  balances 
with  the  energy  absorbed  by  the  molecular  system;  i.e.,  A(a)B(o,T). 
Therefore,  the  emitted  energy  by  the  molecular  system  to  the  radiation 
field  is  given  by  A(o)B(a,T). 


B(a,T) 

Thermal 

Bath 

moleculaj:* 

system 


Absorbed  Energy 
=  B(o,T)A(o) 


[1  -  A(a)]B(a,T) 
Transmitted  Energy 


Figure  5 


The  situation  described  above  is  schematically  shown  in  Figure  5.  A 
thermal  bath  is  a  blackbody  radiation  of  temperature  T.  The  moleciilar 
system  absorbs  the  energy  A(o)B(o,T)  and  transmits  the  energy  [l-A(o) ]B(o,T) . 
The  absorbed  energy  A(o)R(o,T)  is  re-emitted  isotropically  from  the  molecular 
system,  thereby  establishing  the  energy  balance. 

The  discussion  made  above  for  a  single  uniform  layer  can  be  extended 
for  the  case  of  multiple  layers.  If  the  system  is  as  shown  in  Figure  6, 
the  energy  B(T^)A,  emitted  by  the  layer  1 
layer  1 


B(T^)Aj^ 


layer  2 

1  1 

B(T^)A^(1-A2)  '  ' 

layer  N 

B(T2)A2 

1  1 

1  1 

'^2 

1  1 

1  1 

1  1 

1  1 

B(T^)A^(1-A2).. 


B(T2)A2(1-A2)., 


Figure  6 
lU 


.(1-A^) 


i 


becomes  B(T^)A^(A-A2)  .  • .  after  it  goes  througli  the  layers  2  through 

N.  The  emitted  energy  B(T2)A^  by  the  layer  2  yields  R(T2)A^(1-A^) . . . ( l-A^^) 
at  the  output  side  of  the  N-th  layer,  and  so  forth.  Tlic-  totaJ  energy 
available  at  the  output  side  of  the  N-th  layer  i.s  calculated  by  summing 
all  contributions: 

N  N 

E  =  [  E  B(Tjj)Ajj(  ,1  (1-A  )  )  J  +  R(Tj,)A^. 

n=l  j=n+l 


The  computation  requires  the  line  strengths  computed  for  these  temperature;:. 
Since  the  calculations  are  for  temperature  lower  than  300°K,  no  inherent 
difficulties  exist.  The  line  strength  at  temperature  T  is  caLcuJated  from 
the  strength  S(Tg)  at  temperature  T^  by 


S(T)  =  S(T  ) 
s 


-E'VkT 


where  Q^(T)  and  Qp(T)  are  the  vibrational  and  rotational  partition  sum  at 
temperature  T. 

There  are  two  programs  listed  in  the  Appendix,  one  developed  for  tlie 
single  layer  computation  and  another  for  the  muli-i-laycr  (up  to  6)  compu¬ 
tation.  The  line  listing  input  is  TAPE2,  and  other  control  parameters  arc 
in  TAPE5-  The  output  spectrum  produced  on  TAPEl  in  (l6l5)  format  is  a 
normalized  eraissivity  with  respect  to  the  last  layer  (lOOO  for  a  unit 

p 

emissivity).  The  program  SPDG  produces  a  spectrum  smeared  by  sine  .  The 
synthetic  spectra  shown  in  our  scientific  report  were  calculated  using 


these  programs. 


Append i X  A 


Hourcf'  Prop, ram  TREAD 


pK't)( iR AM  r Rt: A 1  T.I  <  r are; V .  tape  i »  f  ape6  » r ape2  » ou  rpur  > 

niMENSION  IA(204U>  ,IC(204{J) 

(tlMENSION  .ll.i(25<i>  .IE<512) 

100  FORMA! (8 ( IX » 08) ) 

101  FORMATdAIS) 

DATA  1XA»  XXB/777777B» 17000000B/ 

N~104 

10  CONTI  NOE 

DO  SO  I==l»256 
ID<  I  )=^0. 

SO  CONI I NOE 

BlIF  FEF<  I N  (  9 » 1 )  (  ID <  1 ) » I D  ( N )  > 

1F(UNIT(V))  llf21»31 

11  CONTINUE 

CALL  TPDCUDEdDflE) 

WRI  I  E  <  A  )  ‘  1:E  ( I )  » I"=l  ,25A  > 

00  TO  10 
;.!J  CUNT ;[  NOE 
;n  CONTINUE 
REWIND  A 
,10=0 
J-=l 

;:0  CONTINUE 
JA='I 

00  15  Is==lv8 
JB==  .)Af2SS 

RF  AD(  A)  (  I:A(  1 )  »I=JA»JB) 

JA==JBF1 
I  S  t:ONTi;NUE 

01.)  12  ;[=•••  I » 20  4  8 

:i:x=  1A(  T  )  .AND.  IXA 

1Y=  ,I:A(  1  )  .AND.IXB 

I F  (TV, NE . 1 5000000B )  GO  TO  12 

:i:c(  j)"ix 

J==Ji:l 

12  CONTINUE 
.)J==J- 1 

U<.J.l.E,l)  GO  TO  20 
II  (  J.GT.  1  )  JJ=--J--1 
JX=  f  <:  JJ/1A))((1A) 

IF  (  IX.GF:  .  lA)  WRITE(2»  101 )  acd  )  I==l  r  JX) 

J0==  .10  + JX 
I  Y~.)X  1 1 
.)  =  I. 

DO  i;!  I  jr»,)j 
i(::<.,))==ic(  I ) 

J  f  1 

IJ  CONTINUE 

II-  (  10 ♦LE- 8 19201  GO  TO  20 
<Ll  CUN  T  I NOF. 

FRINF  101 »  JU 
CALI.,  EXIT 
STOP 
END 
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Ccp7  c.  ^  ■ 
peimit  fjl’y  i 


■p  o 

. .  lion 


SUHROUi  INt:  Tf'I.iCODE;(  ID»IE) 

1.1  ,i:  liENS  .1  uN  I D  <2),  .f  F.  ( 2 )  r  TA  <  6  )  ,  IK  ( 5 )  ,  IC  (  4  ' 

DATA  ( If'K  1 )  »  I  -  1 »  A)/77777777O00O00000O00K»7777777.:^0000Br7777Kr 
:l  777700000000()00000()OB » 7777777700000000K ,  77777777K/ 

DATA  (  IC';  i; )  .  1=1  *4)/17B>1700B»170000Ky  17000000B/ 

I  ^ 

.1=1 

DG  11  1=1 fSJ 
I \  ■•.i:d(i.  >  .and.  ia(  1 ) 
iD<  1  I  r<  iXr.  -ja; 

]:X=lDa.  >  .AND.IA(2) 

JD(2;=SHll  r ' iXy  12) 

1x=.I:D<L  )  .AND.  lAC?) 

:i.  r^IDd.f  1 )  .AND,  IA(4) 
lZ2=SHir T( lYy  -48) 
lZ=SHlf'  T(  IXy  12; 
iD(3)  =  izzf:i:z 
1X= IDd.il) .AND.IA(5) 

1B<4,>=SHIFT(  :i:Xy-24) 

n :<5)  =  :i:d(l+i  »  .and.  ia(A) 

l.=Lf2 

DO  12  K=ly5 
i:X::IB(K),AND,JC(4) 

1Y  =  IB(K)  .AND.  ICCT) 
rZ  =  lB(K'  .AND.  1(:(2) 

.t7.7=ID<  K  )  .AND.  :I:G<  1 ) 

lf;:(  ..D^GH  i:(- 1  Cl  Y,  -12)  ISHIFK  IZy-2)+SHIF  K  IZZyBM  lX 
J=Jf  1 

12  CONTINUE 
11  CONTINUE 

IX  =  IDa...)  .  AND.  ]:a<  1 ) 
iB(i  ):-SH:i:FTn:xy-3A) 

I>:=IB(  1  )  .AND.  .i:c(4) 

I Y=IB( 1 ' . AND, ICC G) 

IZ=ID< 1 ) .AND. IC(2) 

I  /  z=:i  B( :i ; . and.  ic(  1 ) 

1 C :  <  J  )  ■■■■  h  H I F  T  f  1:  Y  .  - 1 2 )  +  S  H I F  T  ( I Z » -  2 )  +  S  H I F  T  ( I Z  Z  r  8 )  + 1 X 

RETURN 

END 
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Appendix  B 


Source  Program  SKXX 


I  KOORAH  IT ( TAPES » OUTPUT »  TAPE2=0UTPUT  f INPUT r  TAPEl ) 
li|;Mt:Nai0NIA<2048) 

.100  F0RfiATa615) 

REAU ( a  » 1 00 ) < I A  < I > » I  =  1 » 2048  > 

J=0 

no  10  I  ==1.2048 
Tr<IA(T).LE.0>  GO  TO  10 

10  CONItNUE 
J.T-204H“J 

,10==^  <  800+JJ)/lA))KlA 
PRINT  lOO.JJ.JO 
JK==(:.U)0».J,l -An 
..il==>jK+i;'3 

UR  I  Iff  (2.  1.00)  (■  IA(  I )  .  I=JK.  JL) 

K'l  All  100.. .IX. Jr 
J2"J(M(  JX- J))tflA+JY 
:i  8IIHA~0 

.1-  2048 
K-49 

III)  1.1  :i:“t.i2B 
.LSUMA  I  SUHA  I  IA(K) 
rsurtp  asurtBf IA(  j) 

Is  -K  i  1. 

.I==,l  t 

11  (DNIfWIP 

.'lA'v'i:  f  r  LiiA  r <  isuma  * ) /128 . 

KAOt::  ( f-  LOAT  (  ISUMB )  >  /128 . 

J  I  I  tx-:  BAOE) 

DU  1.2  i:  4y.j7 
X  •f■l..  OAT(  fACD) 

>o--  '.  X  AA0E))K5.64 
IT'<  1  >  -  I.F  IX(X) 

12  CUN  11  Nt IP 
J.'7==J7M 

III)  1.1  T:  IZZ.2048 
lAi 1 )  =1 A< I ;  lb 
12  CONTINUE 

UR  1  (  E  <  1 . 1 00  ;>  <  I A  (  I )  .  I  -49 . 2048 ) 
iK.i  1.4  K==1.25A 

RPAlii  8  .  i,00>  ( lAl  I ) .  I"1 .2048) 
no  lU  1-1.2048 
I  O'  !.  ;-lA«  f  )-in 
10  CONTINUE 

Wl?  I  f  f-  (  I  .  I  00  X  T  A  ( I  )  .  1  =■  1  .  2048  ) 

14  CUNIINIIE 
CALL  f'XIT 
rop 

i.:Nn 
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Souri'i'  Prof'ram  to  oalculato  a  ;\vnthpl.i(: 
Appendix  C  spept.rum  (a  ninf'le  layer) 


PFiOGRAM  bPSPLDT  '  I  APE.  1  =^200B » TAPE2^^200B »  OUTPU  I  ^^^20()B  v  TAPE:  5=  200P  » 
^^:dHJTPUr»TAPE::8=:200E  • 

:EN(E;:GE;:F^  iik a. >yriii 

l.iinE::N:.ilUN  A(30/2,'  *  I  A(2048)  »£«(  1.024) 

DlMEi-NSlON  D(3)  yDIKSl 
DIMPNGIGN  GCSr<7) 

LilfIPNSlUN  E<0lJi,7> 

E::qu  r  oai..i::ni:;E;:  <  a  <  j.  >  *  i  a  <  i ) )  y  <  a <  204'?  > »  b  <  i.  >  > 

:l.  0  0  F-  U  FC  M  A  I  E  F-  :l  t) .  3  .»  E  .1. 0 . 3  y  F  5 . 4  y  P  1 0 . 3  y  2  A  9  y  3  A  3  y  A  3  ) 

.102  FOFcMA^rU.lb) 

1 0 1  F-  0  FF  H  fl  r  ( I '  I.  ('i .  1 0  y  4  (  1 6 . 4  ) 

:F03  FDFO-tAii.  "SpfctFRUh  IP  CALGLIF.  A  ThB  FROM  ‘  F  :l  0 . 4  ■  "  CM  .1.  i  IJ  “  y 
lF.F0,4y  •  CM'l  y  "  ) 

103  Fl)FCMAT':4F:F0.3) 

100  FURMAY  (  7i: 

DATf,  U.:Bl  .i: )  y  .1.  1  -  7)/3Fl  1.3H  2y3H  Oylll  { y  .>Fl  by. .11  (.y3H  7/ 

llA  I  A  (  F<l..l  I  ,  .l  \  F  '  y  1 "  1  *  /  '  ‘  / 1  •  Ij  y  1  ♦  ^)  y  1  4  Ivi  *  1  4  0  »  1  -  0  >  I  <  3  '■  1  *  e  , 

1415926 
I  T ^^^^230  4 
1NI;I|::..3  -4  ' 

N'- 1 
G  ~  4  3  2 

ALN2.'--AL..UG  \  2  4  .» 

AX---v'.12 

T"2vG4 

RE Al.i  <  5  y  1  06  )  (  OCG  f  <  1  )  y  1  ~  1  y  7  ; 

■)  I  T..y(TT....j  >/(  !  )«  I  I  ) 

P  SUM  ==0.1 

R  E!  yA  D  (  5  y  1  01)  W  X  y  W  y  W  At  y  W  B 
REA»(5yl02)  M 
WX==WX*(-1.) 

WD==WX)F!512, 

U)E==W£t)K2  4 

WDD==  <  ( WB+WX  >  ■  WA  •  WX: )  )  /  ( 2  4  *UE  ) 

JEND=IFIX^WnB)F•l 

WXA==UA--Wli 

UIXB==WXAiWE.)K2. 

WXC-WaBFWD 

AX==WX 

BO  10  l=ly3072 
A(1)==0. 

10  CONFINUE 

DO  20  J==lyJENi:i 
lF\INriEx.F:0.2)  GO  fO  17 
1N1IEX:=:0 

11  GON'ilNUE 

F\'E:  All  E  2  y  I  00  )  F'  y  0  y  <  B  <  1  )  •  1  ==  1  v  21  )  y  E  I B  (  1  )  y  1  1  y  6  > 

BO  A  l==ly7 

6  IF-  E. 111(6)  .EI>4  lii:i(.l>)  l.yri  To  / 
l.iO  I  I.)  I  1 


I 


TAPE6 


/'  1  > 

UCONSi  ^^UCSKl  ) 

II-  <|■■,l...l  .UXA)  &U  (O  11 
IFd'.Bi.WXl)  1X1*1:  X=^i 
If  (N.NE.O)  Wf<ITE(6»  10<)>  l  yLl 
IV  CONTINUE 
QA~Q)(ctJCONS  1 
t!X:^:<  i  .4J9  ))K  ri  r»Ii(2) 

OR  I  (  I  T/1  )  ♦♦A.I 

OA  :  OA)!(PSUM*OR  I  »  t  EXP  <  OX  ) 

OC" o*p;.iUrti«ok  i  * < exp  v ox  >  ,• 

WRiiE(a»ioo>  p).ac*<j;ui;»i-:  j  ,2,  v  '.iiui; 
Ul>  :OA/Wx 

II  '.N.NE.O >  WRHE<6>105)  OyOA»Ol.i 
PIO::(P-^WXAl/UX 
IPR==IT  IX(PR>  i  1 
A( IPR * -Af IPR /  »  OB 
12  CONTINUE 

1I"<1NDEX.NE.;1>  GO  TO  11 
1.J  CUNllNOE 

BO  1!^  [:ay2U'1B 
X-:=A(  I ) 

If  (X.LE.  .  I'J)  GO  TO  14 
If  (X.GE.740.  )X==100, 

IfXX.EE.  -6X0.  }  X=~100. 

A  (  i:  ) 1  .  -  1 1 .  /EXI"- 1  X )  ) 
l-'l  X"::lOOO.>KA(  I ) 

.IA<  1  )  'IF-  1X<X  ) 

JO  CONIINUL 

W  R  I  1  E  (  6  »  1 X  >  Ul  X  A  »  U  X  B 

WR I  (  E  <  1  f  102  )  (  I A  <  I  )  r  I  "Cl  y  2048  ) 

II""  1025 

BO  j6  I-'lyI024 
A  (  I  )  '-"B  ( 1 1 
A  (11/ ""() . 

B  ( I )""(). 

1 1  ""II  f  -l 
16  CONllNUE 
GO  (0  19 
21  EONFINUE 
INDEX  "=2 
00  ro  15 
IV  CONllNUE 
WXA"  N>;  u 
WXBWXAfWEKi:. 

UXl,  ■  I*  /  B  t  WB 
20  CONllNUE 
CALL  EXn 
0  COP 
END 
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Append! x  D 


Rnurcf  I’rofifiiin  ir,  c.-il  I’ulate  a  syntdit^tic 

sped  rum  ( (1-]  ay 'f  ) 


I  'ROGh'AM  MCa:.!Cil:::N  (  TAPE  1  ~20u£<  r  TAf ‘T' :?  -  y  OU  i  I  U  I  .‘OOB  y  TAPE:.;  ■  ■20<)fc'  r  TAf  EA 

f  PU  I  y  TAPE (<==^20010 

THI.3  PKuGKAA  IS  moheeiei.i  i  Rurt  t;u2(H:R  (wh:u;ii  is  GlKin  FOR 
ONE  LAYER)  y  ANIi  USE  FOR  TO  CaILUTAJE  6  l  ArLi:  A/rtUSPHERE 
EACH  Layer  has  HEFFERENI  IEMPERhIORE  aNM  lHliLKEMl  i;uNl..l:  N  I  RA  I  .1  "(V) 
OF  H20y  C02yG3yN20yC0yCH4yO2 
INTEGER  11.1(6)  yii  1.1 

II  :r  ME.NS  1  (.)N  A  <  6  .  1 2  ;  »  a  a  (  S 1 2  >  r  l  a  <  S  l.  2  ) 

1:1:1:  HENS  I  ON  1.1(0  .' ylil.KO) 

DIHENSION  UCS  I  (OyO)  yllB(b')  y  Vl  (0)y  El  I  (0) 

DIHENSION  ROI  JiOly  Tdi) 

EOU I OALENCE  (  AA  (  ;l  >  y  .1 A  (  I  )  ) 

:L  ()  0  F  0  R  M  A  I  (  T  1  (.' «  .1  y  t .  :l.  0  *  J!  y  F-  S  »  4  y  f  1 C’ « .5  y  ••.'  A  V  y  .J  A  8  y  A ) 

102  FORMA  T(  161  IT) 

1 0  :l.  f-  0  R  I'l  A  I  ( F- 1 6 . 1 0  y  •!■  F  :l  A .  4 ) 

10;5  FORHAK  “SPECTRUM  IS  CALCULATEi:'  FROM  “y  FIO,  |y"  Lh  :t  TO“y 
IF  10. 4  y  ■  CH'l  .  ”  ) 

105  FORMAT ( 4E10 . S ) 

106  FORMA  T  (  7t:iO  .  3  ) 

10/  FORMA  T  (  OF- 6 . 1 ) 

DATA  (DTK  1  )  y  l'=l  y  7)/3H  1  y  3H  2y3H  3y3ll  4y.'SFI  Sy.'ll  6y3FI  7/ 

DATA  (ROTJC  I  )  y  I-El  y  7  )  / 1 , 5  y  1 . 0  y  1 . 5  y  1 . 0  y  1 . 0  y  1 . 5.1 . 0/ 

DATA  ClyC2/l  ,  ly()9E-12yl  .4338/ 

PI-/T.  1415726 
INDEX^K) 

N=1 

G=''.32 

AEN2'"A10G(2.  ) 

AX'--0. 12 
TEP  ■:-295- 

READ  (5.107/  (TTd/yt  1.6) 

DO  1  l  -'El.  .6 

R  E  A  Ti  (  5  >'  C'  6  ’  (  Q  C  S  I  *.  f  y  2 )  y  .J  ="1.7  ) 

1  CONTINUE 

READ (5.101)  WX  y  W . WA  y  WB 
READ(5yl02>  N 
DO  3  l-l.o 

3  T  I  T  ( I )  ^=--  (  T  T  (  I  )  --  TEF:-  )  /  (  rE:P* TT  (  I  )  ) 

WR I  I  E  ( 6  y  1 07 )  (  T  r  ( I )  r  I  1.6) 

WR I  I  E  (  6 .106)  (  (  QCS I  ( I .  J  >  .  v)==  1 . 7 ) .  I  ■“  I  y  )> ) 

PSUM==0. 1 
W./==W ;<*(  •  !  .  ) 
wD==wx:xeo. 

WE==IJTK<2. 

WDD==  (  (  WB  F  WX  :>  (  WA  •  WX  )  )  /  ( 2  .  JKWE  > 

JEND==IFIX(WDD)I1 
WXA==WA  - WD 
WXB==WXAFWE)K2. 

WXC- WXBfWD 
AX==WX 

DO  9  J==ly3 
DO  10  I==ly512 
A ( J  y I ) =0 . 
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.1.0  CONTINUI:: 

9  CGNTTfvfUE 

IiD  ('■0  ,)  ”  .l  »  Jt'.No 
INDEX^^^O 
1,1  CONTINUE 

111- 100)  py o»  < 0(  i ) » ]>a  »2)  y  <  iiK  i;  >  y  i>=:i  r6) 

.i  r  U'.le.wxa)  bu  TO  11 
il-  <  F'.OEyWXb)  INDEX: -1 
ilO  A  1  I  y  7 

A  :i:  i-  1 D  .  t  o .  i m  <  :i:  >  >  oo  i  o  / 
bu  10  1.1 

H..,i  ko  1  j<  i: ) 

I  i  I.'  N  ,  Nl::  (  b  I  Ull<  i  i  E  t  6  »  I  00  >  F‘  r  0 
DO  19  KaIvA 

Oi;ONS  1  "OiCb  Ky  I  » 
oA:-:i.T*u(::oNy  i 

OX-"  ^  1 . 4AV  ri  i  I  I  <  K  'XiD  lil!) 

Ok'i':"  (  T  I  (  K  )  /  1 1  I-  .)  .ioEAJ 
OA-OA*kbUfUOR  I  KlEXEbOX  )  ) 
i.D 0  )i<  I-  ■  0  u  f'l  *  0 1  >  i  k  •  e  x  f  •  •.  o  x .» -■ 

1 1'  I\  *  l:,0  •  J  '  OF'!  1  I  i::.  (  0  y  I  v)0  )  F'  y  OL  y  ^  lb  1  *  y  1 I.  y 1  1 D  I. 
OD:-:QA./WX 

!l.  I"  '•  N  i  NE  <  0  )  WIa'  .1  ( F;.  (  6  r  i  (>'::> )  0  y  U A  t  OE* 

F-k-vl'  WXA).--WX 
ir'k-i:i- iA(F'k>  i  i 
A  ^  K  1  1  l"'k  ;  "-I'l k  V  1 1  'k  '  1  Oil 

19  i:;ONi.FNUF- 

II  '>  t  I  IMV.b  I  .  900)  .and.  <  l■•.GT.  WXD)  )  INDEX"-! 

12  CONilNUE 

IF  (  I  NDI::  X  -  NE  .  .1. )  00  ro  11 
.1  s  Cl. IN  I  I  NI..II: 
iC'iAV  ■:  ,1  |■■■k 

i  'K  I.  N  I  X  «  "  I  MAX"-  ”  y  IHiTX 

DO  wb  f:-:  I.  ybCO 

F'-OXh  f-  OX  *  l  i  uAld)  •!  O.SiXWX 
.iH„)  J.  /  I  ,  -  1  V  Cl 
lEriF  "I  I  <10 

.DLACK  DODi'  '.F'LAN(,I\  FUNCTION)  l  UNOH.lN 
ClyC2  are  DOMF;:  CONSFANIy  F-'  Ib  WAOE  NO.  I  EMF'  I  ! 
liD<  k  >:-ClK<  F'K*b  )*!  .0/(EXP(  <C2*P)/TEMP)  --1 .0) 

/  f'l  (  i\  I.  ) 

11  <  A  .  I.E  -  10  )  oo  10  20 

IF  <X  ,ivE.  /'-(().  >>;r:10<), 

11-  <X.LE.  -A/0.  )  X  -  lOO. 

00  10  2 A 

.'b  l(l\,'-1.0  X 

00  lO  2/ 

'.’/i  I  '  i\  )  1 . 0  '  I  .<l''<  X  >  > 


t-'Y 

permit  t  ;i/  i 


>  y  I~PI  y6) 


I  ErtPEKA  I  Ul 
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27  CONTINUE 

AA(l)==(l~r(l)>#£tB(l) 

2  +T< 1 )*( 1-T(2) )*BB(2) 

3  +T(  1  )*T(2)*(l--r(3)  )»BB(3> 

4  +T  <  1 )  *1  C  2 )  ♦  r  (  3 ;  )!i  ( 1  -  r  1 4  >  >  !«BB(  4  ) 

5  +7  ( 1  )  *1  <  2 )  *T  (  3  )  ♦  r  (  4  )  *  ( :l  -  T  (  5  )  )  #B7:f  <  S  ^ 

6  +T  ^  1 )  *T  ( 2 )  *T  ( 3 )  *  T  (  4  >  *1  ( S  )  *  1 1  •••  T  (6 )  )  *riB  (  6  ; 
XX:^(  1000.0  ♦  AA<I))/BB(1> 

IA(  I  >==IFIX(XX) 

35  CONTINUE 

URITE(A,103)  UXArWXB 
WRITE-  ( 1 , 102)  ( IA(  1 .)  f  I=a  »320) 

IL==IMAX  -  320 
DO  45  K  =  l»<fi 
II=-1L  +  1 
JJ^321 

DO  46  1=1 jIL 
A ( K . I) =A ( K » J J ) 

A ( K  f 1 1 ) =0 . 0 
II  =  II-H 
JJ=JJ-H 

46  CONTINUE 

DO  47  I=II»512 
A<K» I )=0.0 

47  CONTINUE 
45  CONTINUE 

UXA=WXB 

UXB=WXA+WE*2. 

UXC=WXB+UD 
60  CONTINUE 
CALL  EXIT 
STOP 
END 
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Appendix  E 


Soui’ce  Program  SPDG 


PROGRAM  SP(3RD  (  TAF’E!  =200ti » TAPE2==200B  f  0UTPUT=200Ei » TAPE6»0UTPUT ) 
DIMENSION  IA( 1024) tA<4096) »8(S12) y IN0<512) »M(3) 

DIMENSION  10(128) 

COMMON  lAySylNV 

100  FORMAT! 1615) 

101  FORMAT! IH  ) 

M!  1  )==11 

M(2)^0 

M!3):=0 

LM-=1 

./512, 

RLAIi!  I  y  100)  (  IC(  I  )  >  I“1  y  128) 

IF  (EOF(l)  .NE.  0)  00  TO  60 
DO  50  L=^lylOO 
DO  9  1=1 y 128 
lA! I )  =  IC! I ) 

9  CONUNUE 

READ ( 1 y 1 00 ) ( I A ! I ) y I  =  1 29  y 1 024 ) 

IF  (EOF! I)  .NE.  0)  GO  TO  60 
J=897 

DO  8  1=1 y 128 
!(.;<  I  )  =  1A!  J) 

J  =  JT  1 

a  CONUNUE 
K=2047 
J--1024 

DO  10  I=Iyl024 
IX=IA( J) 

A!K)=FLOAT( IX) 

A<KH  )=0. 

K  =  R"-2 
J=J-1 

10  CONTINUE 
A!2()49)=0. 

A ( 2050 ) =0 , 

J=4095 


DO  11  I =3 » 2040 » 2  - 
A  ( J )  ==A  ( I ) 

A(J+1)=A(I+1) 

J=J-2 

11  CONTINUE 

CAL  L  HARM  <  A » M » INO . S » 1 » lERR ) 
AX-'=1  .-XX 
J  =  4()95 

DO  12  I=3»1024»2 

A( I)^A( I )»AX 
A( J) =A( J)*AX 
J=^=J-2 
AX-AX-XX 

12  CONTINUE 

DO  13  I “ 1025 » 2048 » 2 
A(  I  r==0. 

A(  J)==0. 

J-J-2 

13  CONTINUE 

CALL  HARM ( A » M . I NO » S » - 1 » lERR ) 

J:=:l 

DO  14  I~l»2048»2 
IA(J)=-IF"IX(A<I)  ) 

J^J+1 

14  CONTINUE 
IFd.NE.l)  L.M-65 

UR I TE ( 2 » 1 00 )< I A ( I ) » I -LM » 960  > 
50  CONTINUE 
60  CALL.  EXIT 
STOP 
END 
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b(?lng  =  1).  A  presence  of  the  N^O  band  In  the  observed  data  Is  rather  questionable 
and  It  will  bo  addressed  when  Irr.proved  data  becomes  available  from  the  future  flight. 

The  measuromont  during  the  Initial  flight  war.  not  carried  out  at  the  optimum  conditions. 
Tb.e  slgnal-to-nol3f'  ratio  was  down-graded  by  unexpected  instrumentation  problem.*;.  Even 
w:ih  the  unfavorable  conditions,  the  data  provided  more  than  ample  evidence  that  the 
experiment  Is  extremely  effective  for  atmospheric  study.  The  project  was  and  will  be 
carried  out  under  a  code  name  of  SCRIBE  (Stratospheric  Cryogenic  Interferometer  Balloon 
Experiment ) . 
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